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Solid Insulating Materials as well as for 
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Department of Electrical Engineering, 
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Abstract 

Interfaces consist a most important part of conventional 
insulating systems at high voltages. They are considered to 
be problem areas which have to be dealt with. Numerous pub- 
lications have contributed in rendering the mechanisms of 
interfaces understandable. On the other hand, interfaces in 
nanocomposite polymers seem to function in an entirely dif- 
ferent manner from that in conventional insulating systems. 
The present paper reviews past work on both the conventional 
insulations and in nanocomposites. Differences regarding the 
interfaces are mentioned and discussed. Whereas interfaces 
in conventional insulating systems are to be avoided, inter- 
faces in nanocomposite polymers seem to be desirable - at 
least - up to a certain percentage of nanoparticles in the 
base polymer. Although things are better understood in con- 
ventional insulating materials, more work has to be per- 
formed in order to clarify several aspects, such as space 
charges and electrical trees emanating from enclosed cavi- 
ties. Needless to say that much more work has to be done in 
nanocomposites w.r.t. their modeling and possible explana- 
tions of the surprising performance of interfaces, a perfor- 
mance that deviates strongly from the performance in classi- 
cal insulating materials. 



Keywords 

Breakdown, breakdown strength, conventional insulating 
materials/systems, nanocomposite polymers, partial dischar- 
ges, electrical treeing/trees 
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Introduction 

With some sort of exag- 
geration, a well known pro- 
fessor said once that "the 
problems of high voltage in- 
sulations are problems of in- 
terfaces" [1] . Interfaces re- 
sult when there are two dif- 
ferent insulating materials 
next to each other or when an 
insulating material meets a 
conductor. An interface may 
become, e.g., the source of 
partial discharges or even the 
cause for a complete break- 
down of an insulating system. 
When two insulating materials 
of different dielectric con- 
stants have a common inter- 
face, then the material with 
the lowest dielectric con- 
stant will undergo the more 
intense electric stressing 
[2] . Depending also on other 
parameters, such stressing 
may result in a gradual dete- 
rioration of the insulating 
system and consequently in a 
complete breakdown [3] . 

In this paper, the ques- 
tions raised by the existence 
of interfaces in high-voltage 
insulating systems are dis- 
cussed. Interfaces play an 
important role in determining 
the robustness of an insulat- 
ing system, when conventional 
insulating materials are 
used. Interfaces play also a 
vital role in insulating sys- 
tems with nanocomposite poly- 
mers. Their functioning, how- 
ever, is of another nature. 



This short review is by no 
means exhaustive since the 
topic of interfaces is a vast 
one. It is the aim of the au- 
thors to give the gist of the 
problems and questions the 
researchers may face and to 
offer some comments. In the 
context of this paper, the 
terms "insulating material" 
and "dielectric material" or 
simply "dielectric" are some- 
times interchanged meaning 
the same thing. 

Interfaces in Conventional 

Insulating Materials 

Dielectric breakdown in 
insulating materials depends 
on electrode configuration, 
insulating material thick- 
ness, electrode materials, 
presence of microcavities, 
temperature, pressure, nature 
and morphology of material 
under test, type of applied 
voltage, damage path (surface 
or volume) [4] . Various di- 
electric breakdown theories 
have been put forward [5]- 
[13] . No matter whether the 
proposed theory was based on 
cumulative impact ionization 
by electrons - creating thus 
positive space charges which 
distort the field distribu- 
tion and weaken the dielec- 
tric [5] -, on the notion of 
"intrinsic breakdown" - ac- 
cording to which a large num- 
ber of electrons trapped in 
energy levels due to lattice 
imperfections can transfer 
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energy to the lattice vibra- 
tions [6] - [8] , on the "40 ge- 
nerations avalanche theory" 

[9] , on the importance of spa- 
ce charges which modify the 
local electric field value 

[10] , on the right assumption 
that the breakdown is a prop- 
erty of the dielectric mate- 
rial plus its electrode sys- 
tem [11] or on the theory 
based on the ionizing elec- 
trons and the hole traps 
[12], [13], the fact remains 
that all the above mentioned 
phenomena result from elec- 
tric field intensifications, 
i.e. from either electrode 
imperfections or mismatch of 
dielectric materials. This 
brings us to the point men- 
tioned in the Introduction of 
the present paper: that in- 
terfaces may create the con- 
ditions which may cause elec- 
tric field intensifications. 

The subject of partial 
discharges (PD) which may en- 
sue because of electric field 
intensifications and/or be- 
cause of gas (or foreign par- 
ticle) inclusions in a solid 
dielectric material, has been 
studied by Mason in his fun- 
damental publications [14]- 
[18] . Having in mind all the 
above, it is fitting to say 
that interfaces - created ei- 
ther by a mismatch of dielec- 
tric materials or because of 
intrusion of foreign parti- 
cles and/or air cavities in 
the insulating material under 



question -, are the problem 
areas of an insulating sys- 
tem . 

As was pointed out quite 
early [19], interfaces play a 
most significant role in the 
discharge and breakdown pro- 
cesses: even if an insulation 
does not contain any cavi- 
ties, at a sufficient stress 
"some event" releases gas 
with the subsequent formation 
of a cavity. The cavity is 
occupied by a gas discharge 
which increases the rate of 
gas formation with the subse- 
quent growth of both the cav- 
ity size and the discharge 
intensity. The importance of 
the differing nature of in- 
terfaces was also stressed in 
another publication, where it 
was pointed out that damage 
in internal cavities in poly- 
ethylene is little compared 
to the electrode adjacent ca- 
vities of the same dimensions 
and tested under the same ex- 
perimental conditions [20] . 
It is evident that in [20], 
interfaces between polyethy- 
lene and gas were compared 
with interfaces between poly- 
ethylene and electrodes, and 
the latter were found to be 
more dangerous and deleteri- 
ous to the insulating mate- 
rial. On the same lines, 
Kreuger showed that with PVC- 
insulated cables, the number 
of discharges increased with 
increasing electric stress in 
the dielectric [21] . In yet 
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another paper, it was indi- 
cated that the nature of in- 
ternal discharges was greatly 
affected by the assembly of 
the electrode system and the 
adhesion of the insulating 
tapes [22] . Discharges always 
start in the electrically 
weaker insulating material, 
as was commented in previous 
works [23] - [25] . 

Needless to say that phe- 
nomena related to PD, such as 
electrical treeing, are also 
closely connected to the mis- 
match of the dielectric con- 
stants of insulating materi- 
als and/or to the existence 
of gas cavities in their vo- 
lume. Earlier papers indi- 
cated that the treeing phe- 
nomenon in polyethylene ca- 
bles started from both inner 
and outer surfaces and also 
from solid particles and fi- 
bres [26] . Pioneering work 
performed with 15 kV and 22 
kV polyethylene insulated ca- 
bles reported that trees ori- 
ginated from contaminants and 
cavities, the tendency for 
tree initiation from a con- 
taminant being probably more 
affected by the contaminant 
material than by the size, 
location or shape of the con- 
taminant [27] . The importance 
of enclosed cavities in the 
initiation of trees was also 
reported more recently [28] . 
According to Ieda [29], tree 
propagation can be induced by 
internal gas discharge in the 



tree. It is to be bore in 
mind that in numerous publi- 
cations dealing with experi- 
mental work, the electrode 
arrangement that was used was 
a needle-plane electrode ar- 
rangement, indicating again 
that an electrode arrangement 
was chosen, with pronounced 
interfaces, in order to study 
the treeing phenomenon [30]- 
[33], 

Interfaces which may play 
a role in determining the 
breakdown strength of an in- 
sulating material need not be 
only interfaces between insu- 
lating material and metal or 
between insulating material 
and gas cavity or contami- 
nant. Differing phases may 
play also a role, as was no- 
ted in [34], where the inter- 
facial domain of crystalline 
and amorphous phases may de- 
termine the various proper- 
ties of semi-crystalline po- 
lymers, such as biaxially ori- 
entated polypropylene (BOPP). 
Stressing the importance of 
interfaces and experimenting 
with cross-linked polyethy- 
lene (XLPE) , McKean showed 
that a considerable improve- 
ment in cable breakdown can 
be achieved by impregnation 
with silicone oil or diethy- 
leneglycol. Such liquids can 
impregnate gas microspaces in 
the main insulation and thus 
result in an increase of the 
breakdown strength [35] . Si- 
milar observations were also 
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reported with polypropylene 
and polyethylene impregnated 
with suitable dipolar liquids 
[36]. 

Interfacial breakdown was 
studied with various elec- 
trode geometries and insulat- 
ing systems consisting of pa- 
per typical for transformers 
and transformer oil [37] . It 
was reported that interfacial 
breakdown will occur if the 
paper is not carefully dried 
or if many gaseous micro- 
porosities are left in or on 
the paper. In [37], however, 
it was also noticed that us- 
ing a carefully prepared pa- 
per-oil interface structure, 
the breakdown does not neces- 
sarily take place at the in- 
terface. Similar observations 
were made more recently by 
using silicone rubber inter- 
faces, where both perpendicu- 
lar and parallel to the ap- 
plied electric field were in- 
vestigated [38] . 

Conventional paper-oil ca- 
ble insulation was studied 
quite early and the problems 
of interfaces were noted 
[39] . Alternative insulating 
systems, based mainly on po- 
lymeric materials, were pro- 
posed with considerable com- 
mercial success [40] - [43] . 
Modern cables with solid 
polymeric insulation did not 
avoid the problems of inter- 
faces, namely those of extru- 
sions of semi-conducting 
sheaths with the main insula- 



tion or the inclusion of mi- 
crocavities and/or impurities 
[44], [45]. Extruded cable in- 
sulation exposed to wet con- 
ditions suffered from elec- 
trochemical treeing and impu- 
rities greatly deteriorated 
its electrical performance 
[46] . Moreover, operating elec- 
trical stresses may also cau- 
se premature insulation fai- 
lure in 15 kV polyethylene 
cables, if combined with un- 
favorable interface profiles 
and moisture [47] . Interfaces 
between polyethylene and 
small contaminants or micro- 
cavities may cause bow-tie 
trees in polyethylene cables 
[48]. 

Interfaces either perpen- 
dicular to the applied elec- 
tric field or parallel to it 
or at an angle with it were 
dealt with in [49], where it 
was noted that such a variety 
of interfaces may be encoun- 
tered in applications, such 
as capacitors, cables and in 
transformer windings. Compo- 
site insulating systems must 
preserve low dielectric los- 
ses. Higher dielectric losses 
may imply high ionic concen- 
tration in a solid/liquid in- 
sulating system, i.e. high 
ionic concentration in the 
liquid component of such a 
system [50] . 

Before concluding this 
section, it is fitting to men- 
tion the composite insulating 
systems of electrical machi- 
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nes, which consist mainly of 
epoxy resin and mica sheets. 
Previous work done in this 
context indicated that 
electrical treeing propagates 
through the epoxy resin and 
generally stops at the mica 
sheets, as mica is harder and 
electrically stronger than 
epoxy resin. The importance 
of such interfaces was re- 
ported before using a needle- 
plane electrode arrangement 
[51], where experiments were 
carried out without and with 
a mica sheet inserted in 
epoxy resin (Figs. 1 and 2). 
Evidently electrical trees 
were propagating more easily 
in the case of absence of the 
mica sheet and with much more 
difficulty with mica sheet. 

Simulation work done re- 
cently showed that mica she- 
ets prevent electrical trees 
from reaching the opposite 
electrode [52] . The purpose 
of mentioning experimental 
results regarding the compo- 
site system of epoxy resin/ 
mica sheets is to show that 
the electrical trees propa- 
gate through the electrically 
weaker insulating medium. The 
simulation results indicate 
that even the slightest vari- 
ations of dielectric constant 
may cause the electrical tree 
growth. In other words, the 
simulation data indicate that 



local fluctuations of dielec- 
tric constant imply local - 
even microscopically minute - 
formations of interfaces, 
which in turn may mean local 
field intensifications, en- 
couraging thus the growth of 
electrical trees (Figs. 3 and 
4). Such observations w.r.t. 
the local variations of die- 
lectric constant have also 
been reported for polyethy- 
lene [53] , [54] . 

It is evident from all the 
above that interfaces in 
classical insulating systems 
seem to cause problems (pos- 
sible dielectric constant 
mismatch, PD, treeing pheno- 
mena and ultimately risk of 
ultimate insulation failure). 
Due attention should be paid 
in choosing the insulating 
materials for specific appli- 
cations and to the construc- 
tion of the composite insu- 
lating system. Too many 
things depend on the quality 
of the construction of the 
interfaces [55]-[58], too ma- 
ny things that cannot be ig- 
nored. Interfaces in tradi- 
tional insulating systems are 
considered as the weak as- 
pects of such systems. Keep- 
ing this in mind and without 
exaggerating, it is not far 
from the truth if we state 
that an insulating system is 
as good as its weakest inter- 
faces . 
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Fig. 1: Electrical tree propagation without the presence of 
mica sheet (applied voltage 28 kVrms, 50 Hz) (after [51]) 




Fig. 2: Electrical tree propagation with the presence of 
mica sheet. The mica sheet increases the propagation time of 
the tree (applied voltage 28 kVrms, 50 Hz) (after [51]) 
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Fig. 3: Simulated electrical tree propagation with one mica 
sheet. Needle-plane electrode arrangement used 




Fig. 4: Further expansion of electrical trees. 
Electrical trees stop at mica sheet. 
Needle-plane electrode arrangement used 
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It goes without say that 
this short review regarding 
classical interfaces does not 
by any means cover the whole 
subject and variety of solid 
insulating materials (for ex- 
ample, no mention in this pa- 
per was offered about the in- 
terfaces in outdoor polymeric 
insulators [59], [60] or in 
indoor polymeric insulators 
[61]). Both in the libraries 
and in the Internet, the in- 
terested reader may find pra- 
ctically tens of thousands of 
publications referring to the 
questions and problem areas 
of the solid insulating mate- 
rials and insulating systems. 
What this short review tried 
to do is to show that inter- 
faces, electric field inten- 
sifications, pre - breakdown 
phenomena (such as PD and 
electrical trees) and break- 
down mechanisms are all in- 
terwoven and interrelated. 
Having said that, the next 
question related to inter- 
faces, is whether they play 
the same detrimental role in 
the new generation of insu- 
lating materials, the nano- 
composite polymers. This is 
to be examined in the follow- 
ing section. 

Interfaces in Nanocomposite 

Polymers 

More than twenty years 
ago, nanocomposite polymers 
came to our lives [62] . The 
first application of nanocom- 



posites appeared in 1990, 
when Toyota Motor Corporation 
introduced nanocomposite ny- 
lon in their car industry 
[63] . Since that year, many 
car industries introduced the 
use of nanocomposite poly- 
mers. Use of nanocomposite po- 
lymers was noted in other in- 
dustries, such as in the op- 
tics and the electronics in- 
dustries as well as in the 
food industry. A seminal pa- 
per by Lewis gave the impetus 
for research also in the in- 
sulation branch [64] . 

For the electrical insula- 
tion, nanocomposite polymers 
are defined as conventional 
polymers in which particles 
smaller than 100 nm are added 
and dispersed in such a way 
that at the end one gets a 
homogeneous mixture [63] . The 
addition of such nanoparti- 
cles (the term "nanofillers" 
is also widely used) is being 
done in very small quanti- 
ties, usually less than 10 wt%. 
Nanocomposite polymers con- 
sist of three components: 

a) the base polymer (or poly- 
mer matrix), 

b) the nanoparticles (or na- 
nofillers) and 

c) the interaction zone (or 
interface zone) between the 
matrix and the nanoparticles 
[63], 

Regarding the polymers used, 
these may be either ther- 
moplastics, thermosettings or 
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elastomers. Nanoparticles may 
be classified w.r.t. their 
dimensions, and they can be 
distinguished as 

1) mono-dimensional (i.e. ex- 
tremely thin), 

2) two-dimensional (nanotu- 
bes) and 

3) three-dimensional (inorga- 
nic oxides) . 

The most usual nanoparticles 
for the purposes of electri- 
cal insulation are 

1) silica nanoparticles SiC> 2 , 

2) montmorillonite nanopar- 
ticles (layered silica), 

3) metallic oxides such as 
AI 2 O 3 , Ti0 2 , MgO and ZnO and 

4) carbon nanotubes. 

Nanocomposite polymers can be 
obtained in two types of 
structures, namely, 

(i) intercalated nanocompo- 
sites (formed when there is 
limited inclusion of polymer 
chain between the clay layers 
with a corresponding small 
increase in the interlayer 
spacing of a few nanometers 
and 

(ii) exfoliated nanocomposi- 
tes (formed when the clay 
layers are well separated from 
one another and individually 
dispersed in the continuous 
polymer matrix [65], [66]. 

As mentioned above, nano- 
particles are added and dis- 
persed in relatively small 



quantities in the base poly- 
mer (usually no more than 10 
wt%) . Since nanoparticles are 
smaller than microparticles 
(smaller by three orders of 
magnitude), their interaction 
with the surrounding polymer 
matrix is much greater [62] . 
The so-called interaction zo- 
ne is the main factor con- 
tributing in the improvement 
of the properties of the base 
polymer [67] . In the case of 
addition of nanoparticles in- 
to a polymer, the interfaces 
are far more numerous and far 
larger than in the case of 
microparticles. As the size 
of the added particles is re- 
duced, the interface becomes 
larger and larger. The dis- 
tance between the nanoparti- 
cles is also extremely small. 
It seems that interfaces de- 
termine to a great extent the 
properties of nanocomposite 
polymers . 

The size of nanoparticles 
and the distance between them 
is of the order of magnitude 
on nanometers. Such particles 
may interact with the polymer 
matrix both physically and 
chemically in the nanometer 
scale. This has as conse- 
quence the appearance of pro- 
perties that are somehow dif- 
ferent from those we already 
know in a more macroscopic 
scale [68] . In contradistinc- 
tion to the interfaces in 
classical insulating materi- 
als, and also to what we know 



FUNKTECHNIKPLUS # JOURNAL 



Vl—16 



ISSUE 4 - YEAR 1 




INTERFACES IN HIGH VOLTAGE ENGINEERING: A MOST IMPORTANT QUESTION 



from classical high voltage 
textbooks, the improved insu- 
lating properties of the na- 
nocomposite polymers are due 
to 

a) the large surface area of 
nanoparticles, which creates 
a large interaction zone, 

b) the changes in polymer 
morphology because of the 
large interaction zone, 

c) the changes in the space 
charge distribution and 

d) a dispersion mechanism 
[69]- 

Both the size of the nanopar- 
ticles and the chemical pro- 
perties of their surface play 
an important role in deter- 
mining the electrical, ther- 
mal and mechanical properties 
of nanocomposites. Needless 
to say that the chemical com- 
patibility between the intro- 
duced nanoparticles and the 
polymer matrix is of para- 
mount importance for the ge- 
neral properties of the nano- 
composite [70] . 

One of the most signifi- 
cant characteristics of nano- 
composite polymers is the in- 
crease of their breakdown 
strength as the size of the 
added nanoparticles tends to 
extremely small values. This 
increase is not in agreement 
with the conventional wisdom, 
which suggests that as the 
number of interfaces in- 
creases, the breakdown stren- 
gth decreases dramatically 



[2], [3]. Nanocomposite poly- 
mers seem not to agree with 
what we already know for 
classical insulating materi- 
als or systems [71] . Diffe- 
rences in breakdown strength 
between conventional epoxy 
resin and epoxy resin with 
nanoparticles was reported in 
[72] . Such observations were 
also noted later, when six 
different materials based on 
epoxy resin with various with 
and/or micro- and nanoparti- 
cles of alumina/silica were 
tested. It was shown that 
epoxy resin with nanosilica 
particles was the most suit- 
able to obtain high values of 
breakdown strength [73] . 

Addition of the percentage 
of nanoparticles to epoxy 
resin up to a certain level 
favors the increase of break- 
down strength both with a.c. 
and d.c. voltage, as was no- 
ticed in [74] . Why nanoparti- 
cles act in such a favorable 
way, despite the numerous in- 
terfaces they create? The in- 
crease of the breakdown 
strength may be due to 

(i) the increase of the sur- 
face area of the interfaces, 
which somehow alters the be- 
havior of the polymer, 

(ii) the changes of space 
charge distribution inside 
the insulation structure, 

(iii) the dispersion mecha- 
nism, and 

(iv) the changing properties 
of the insulating material, 
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more specifically its volume 
resistivity, its tan5 and its 
dielectric constant. 

It is possible that the elec- 
trons moving in such a nano- 
composite polymer, loose 
their kinetic energy because 
of the nanoparticles. Since 
the distances between the na- 
noparticles are extremely 
small, the electrons cannot 
acquire enough speed so that 
they can contribute to the 
breakdown process. Consequen- 
tly, epoxy resin with nano- 
particles presents a higher 
breakdown strength than con- 
ventional epoxy resin [74] . 

Normally the introduction 
of particles in polymeric ma- 
terials has as result the in- 
troduction of defects and 
subsequently the worsening of 
its electrical properties. 
Nanocomposite polymers seem 
not to obey the above rule, 
as the mechanisms of conduc- 
tivity during the breakdown 
process are influenced from 
the applied electric field, 
the dielectric constant of 
the nanoparticles and their 
number. The combined effect 
of these factors is difficult 
to fully understand at this 
stage and we need more work 
[71] . Similar results were 
obtained with epoxy resin 
with nanoparticles of Ti 02 as 
it presented a much higher 
breakdown strength than con- 
ventional epoxy resin [75] . 



On the other hand, elec- 
trical treeing propagation 
was found to be easier in 
conventional polymers than in 
nanocomposite polymers [76], 

[77] . Even a small wt% addi- 
tion of nanoparticles affects 
in a positive way the elec- 
trical treeing resistance of 
the nanocomposite polymer. It 
seems that electrical tree 
propagation paths go through 
the base polymer and around 
the nanoparticles (experimen- 
tal evidence for this was 
presented in SEM photographs 
published in [76]). Conse- 
quently, the more the nano- 
particles in a polymer, the 
more difficult the formation 
of treeing paths. It seems 
that nanoparticles act as ex- 
tremely small barriers, thus 
preventing the easy growth of 
electrical trees. Electrical 
trees propagate through the 
base polymer (in other words 
through the polymer matrix) 
and not through the nanopar- 
ticles. In some cases the 
electrical trees stop at the 
nanoparticles and they do not 
progress any further. Such 
observations were made in 
simulation studies recently 

[78] - [80] . 

Further research showed 
that a small percentage in- 
troduction of nanoparticles 
into a conventional polymer 
may increase its resistance 
to electrical treeing [81] . 
It is interesting to note 
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that nanoparticles may func- 
tion as barriers preventing 
the tree growth even in 
minute quantities [82] . Ear- 
lier work indicated that as 
soon as the electrical tree 
tou-ches the nanoparticle, 
the physico-chemical proper- 
ties of it are such that very 
high energies are required in 
order to cause its deteriora- 
tion [83]. Although the lat- 
ter paper is an old one, it 
may give a clue as to why 
nanoparticles act as elemen- 
tary barriers and why they 
prevent (or they delay) elec- 
trical treeing. More re- 
cently, a similar argument 
was given by some Japanese 
researchers [84] . 

Loading (i.e. the percent- 
age of included nanoparticles 
expressed in wt%) plays also 
a vital role in determining 
the resistance to electrical 
treeing. More loading (i.e. 
more nanoparticles, that is 
more interfaces) implies a 
better resistance to treeing. 
This is probably because 
trees interact with many more 
nanoparticles and this delays 
their growth [85] . Another 
possible explanation was of- 
fered some years ago, where 
the authors proposed that in 
front of a tree a damage 
process zone is formed in a 
conventional polymer. Such a 
zone cannot progress easily 
when it meets nanoparticles 
[86] . Simulation data indi- 



cated that loading affects 
the tree growth. The nanopar- 
ticle size plays also a role 
in delaying tree propagation, 
smaller nanoparticles offer a 
better tree resistance than 
the larger ones [79] . In 
Figs. 5-8 simulation results 
regarding the loading of na- 
noparticles as well as the 
size of nanoparticles are 
shown. It is evident that mo- 
re loading makes tree growth 
more difficult. It goes also 
without say that smaller 
nanoparticles offer a better 
resistance to tree propaga- 
tion . 

There is no need to empha- 
size that there is also in 
the field of nanocomposite 
polymers a vast body of tech- 
nical literature, too vast to 
be mentioned here in the con- 
text of this paper. From this 
short review it is obvious 
that interfaces in nanocom- 
posites play an entirely dif- 
ferent role from the one they 
play in classical insulating 
materials. Why is this so? 
This may be because the phy- 
sics and/or chemistry somehow 
change in the nanoscopic 
world of such materials. In- 
terfaces become highly desir- 
able - at least up to a cer- 
tain percentage of added na- 
noparticles. The surface area 
of the nanoparticles is huge 
if compared with that of mi- 
croparticles for other high 
voltage applications. 
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Fig. 5: Simulation with epoxy resin filled with HO 2 nano- 
particles (loading of 2 wt%, nanoparticle diameter 100 nm) 




Fig. 6: Simulation with epoxy resin filled with Ti02 nano- 
particles (loading of 6 wt%, nanoparticle diameter 100 nm) 
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Fig. 7: Simulation with epoxy resin filled with TiC >2 nano- 
particles (loading of 6 wt%, nanoparticle diameter 200 nm) 




Fig. 8: Simulation with epoxy resin filled with Ti02 nano- 
particles (loading of 6 wt%, nanoparticle diameter 100 nm) 
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Models regarding the ex- 
planation of the functioning 
of nanocomposites were pro- 
posed as well as how the 
nanoparticles behave inside 
the base polymeric material 
[62], [87], [88]. Such models 
try to explain the higher 
breakdown strength values of 
nanocomposites and also their 
higher resistance to electri- 
cal treeing. The explanations 
seem to be plausible but more 
hard evidence is needed. Such 
evidence will be offered by 
many more photographs (SEM 
and TEM photos) showing in 
detail how the electrical 
trees circumvent the nanopar- 
ticles and how they propagate 
through the polymer matrix. 
Such photographs are very 
difficult to obtain. 

Where all this leaves us? 
How can we understand in a 
unified way interfaces in 
both classical and nanocom- 
posite insulating systems? 
Can the experimental results 
and simulations on electrical 
treeing with typical machine 
insulation ([52]) provide a 
hint also for a possible ex- 
planation in the nanoscopic 
world? A recent paper on in- 
terfaces posed some pertinent 
questions regarding classical 
and nanocomposite insulating 
materials [89] . From this pa- 
per it is obvious that al- 
though more questions are in 
need of an answer for the na- 
nocomposites, the subject of 



interfaces in classical insu- 
lating systems is by no means 
finished. For example, charg- 
ing of larger interfaces, 
such as found in cable joints 
and terminations, needs to be 
further explored regarding 
the mechanisms of space char- 
ges. Another aspect in need 
of further discussion is whe- 
ther electrical trees may em- 
anate from enclosed cavities 
in conventional polymers. The 
latter question has been par- 
tly answered in [28], [90] - 
where some experimental evi- 
dence was offered as to the 
possibility of electrical 
trees stemming from enclosed 
cavities - but further expe- 
rimental data is needed. 

One thing that should not 
be forgotten - and it is com- 
mon to conventional as well 
as nanocomposite polymers - 
is that an insulating system 
to a great extent is as good 
as its interfaces. This means 
that, no matter whether we 
deal with conventional insu- 
lating materials or nanocom- 
posites, preparation and con- 
struction in both cases has 
to be not only careful but 
meticulous . 

A last remark on the li- 
terature presented here: the 
interested reader may find 
that the authors dwell per- 
haps too much in the older 
scientific literature. This 
is not done because they tend 
to ignore the more recent re- 
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search: they simply would 
like to show that even in the 
old days, the problems were 
more or less the same. More- 
over, they would also like to 
show that fundamental ideas - 
which are with us even today 
- on the various mechanisms 
in insulating materials came 
about quite early. 

Conclusion 

This review - by no means 
exhaustive - tackled the sub- 
ject of the importance of in- 
terfaces both in conventional 
insulating materials and in 
nanocomposite polymers. Whe- 
reas interfaces are to be 
avoided in conventional mate- 
rials, they seem to be a 
blessing in nanocomposites. 



Whereas in conventional mate- 
rials they cause problems of 
compatibility and sometimes 
high field intensifications 
with all the bad consequences 
such intensifications entail 
(i.e. PD, trees), in nanocom- 
posites they seem - up to a 
certain loading - to be de- 
sirable and they prevent (or 
the delay) tree growth. Whe- 
reas in conventional insulat- 
ing systems the introduction 
of more interfaces seems to 
cause sometimes insurmount- 
able problems, the introduc- 
tion of interfaces (because 
of the nanoparticles' intro- 
duction to the polymer ma- 
trix) seems to alleviate elec- 
trical trees and to distrib- 
ute more evenly electric 
fields and space charges. 
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Abstract 

The Beltrami fields are extremely important for the desc- 
ription of time-harmonic electromagnetic fields in chiral 
media. This paper introduces the spherical Beltrami fields 
in chiral media, extending the corresponding known results 
for plane electromagnetic waves. Two key theorems are given 
for the scattering of spherical electromagnetic waves in 
chiral environment from a perfect conductor: the Reciprocity 
and the General Scattering theorems. 

Keywords 
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Introduction 

The name "Beltrami field" 
comes from the fact that 
these fields A satisfy the 
equation type V x A = crA , 
where a a 0 is a reciprocal 
length [1] . The word chiral 
regards the concept of "chi- 
rality", derived from the 
Greek word hand. In geometry, 
a figure is chiral if it can- 
not be mapped to its mirror 
image by rotations and tran- 
slations alone. Typical exam- 
ples are human hands, snail 
shells, spirals and spirals 
in general. 

The discoveries of natural 
optical activity in special 



materials have been known 
since the beginning of last 
century. Though optical acti- 
vity has been considered in 
optics and in quantum mecha- 
nics for many years, its 
analysis within the framework 
of classical electromagnetic 
field theory arose much la- 
ter. Recently, there has been 
a considerable interest in 
the study of scattering and 
diffraction by chiral medium. 
Chiral media are isotropic 
birefringent substances that 
respond to either electric or 
magnetic excitation with both 
electric and magnetic polari- 
zations. Such media have been 
known since the end of the 
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nineteenth century (e.g. the 
study of chirality by Pas- 
teur) and find a wide range 
of applications from many 
sciences [1] . 

In general, the electro- 
magnetic fields inside the 
chiral medium are governed by 
Maxwell equations together 
with Drude-Born-Fedorov equa- 
tions in which the electric 
and magnetic fields are coup- 
led [1] . The chiral medium is 
characterized by the electric 
permittivity e, magnetic per- 
meability p and chirality 

measure (3 . 

Scattering problems with 
incident waves have been stu- 
died: for spherical acoustic 
waves in the works [2] -[4], 
as well in the works [2], [5] 
for spherical electromagne- 
tic waves in non-chiral elec- 
tromagnetic waves and finally 
for spherical electromagne- 
tic waves in chiral media in 
the works [6] , [7] . 

In this work, we derive 
reciprocity relation and gen- 
eral scattering theorem, from 
spherical electromagnetic wa- 
ves emanating from point sour- 
ces in chiral media and scat- 
tered by a perfect conductor. 
Similar theorems are proved, 
analogously, for a chiral di- 
electric [6] , [7] . 

As it is well known [1], in 
a homogeneous isotropic chi- 
ral medium the electromagne- 
tic field is composed of left 
-circularly polarized (LCP) 



and right -circularly polari- 
zed ( RCP ) components which 
are propagated independently 
and with different phase 
speeds. When either a LCP or 
a RCP or a linear combination 
of LCP and RCP electromagnet - 
tic waves is incident upon a 
chiral scatterer then the 
scattered field is composed 
of both LCP and RCP compo- 
nents. So, using Bohren de- 
composition [1], [6], [7], 
the electromagnetic waves are 
expressed in terms of LCP and 
RCP Beltrami fields. The use 
of the Beltrami fields for 
problems of scattering in 
chiral material leads to a 
more simplified relationship. 

This is due to the fact 
that the differential equa- 
tions which satisfy Beltrami 
fields is of first order whi- 
le the electric and magnetic 
field satisfy the Helmholtz 
modified equation which is of 
second-order. Moreover, the 
Beltrami equations and the 
conditions of radiation for 
the Beltrami fields are in 
effect separately for each 
one of them but on the scat- 
terer surface, when applying 
the boundary conditions, both 
Beltrami fields are present. 
This fact makes easier to use 
these equations mainly in 
scattering problems where the 
behavior of the scattered 
field away from the scatterer 
is studied. 

In the second section, we 
formulate the problem for 
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spherical electromagnetic wa- 
ve in a chiral medium. In the 
third section we formulate 
the problem of electromagnet - 
tic wave scattering from a 
perfect conductor as a fun- 
ction of the Beltrami fields 
using the Bohren transforma- 
tion. The spherical Beltrami 
fields are defined in the 
fourth section. Finally, Re- 
ciprocity relation and the 
General Scattering theorem are 
given in fifth section. 

Problem Formulation 

Let Q be a bounded and 
closed subset of M 3 having a 
C 2 - boundary S, i.e. dQ=S. 

The set £T will be referred 
to as the scatterer. The ex- 
terior Q + = R 3 /Q~ of the scat- 
terer is an infinite isotro- 
pic homogeneous chiral medium 
with electric permittivity e, 
magnetic permeability p and 
chirality measure (3. 

The scatterer is filled 
with a isotropic chiral me- 
dium with corresponding phy- 
sical parameters £“, p and 
. All the physical parame- 
ters are assumed to be real 
positive constants. 

Let (E r lnc , H* nc ) be a time- 
' '0 ' '0 ' 

Harmonic spherical electro- 
magnetic wave, due a point 
source located at a point 
with position vector r 0 with 



respect to an origin 0 in the 
vicinity of Q . 

This wave is incident upon 
the scatterer Q and let 

(E^, Hj: c ) be the corresponding 

scattered field. Then the to- 
tal electromagnetic field 
(E, 0 , K e ) in CT is given by 

E^(r) = Ef(r) + E-(r) 
f (!) 
\\ r (r) = H“ c (r) + H* c (r) 

r 0 \/ r 0 \/ r 0 \/ 

We assume that the scat- 
tered field satisfies the 

Silver-Muller radiation con- 
dition 

iE-(r) + rxH-(r) = o(r- 1 ), ^ 

r ^ oo 

uniformly in all directions, 
reS 2 , where S 2 is the unit 

sphere in K 3 , r = |r| , r = r/r 
and rj = (p/e) 1/2 is the intrinsic 
impedance of the chiral me- 
dium in Q + . It is known [8], 
that (2) can be replaced by 

rxE-(r)- n H; c (r) = o(r- 1 ), 

r ^ oo 

uniformly in all directions, 
r e S 2 . The total electric 

field E tot satisfies the boun- 
dary condition: 

n x E^r) = 0, r e S (4) 
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where n is the unit vector 
perpendicular to the outer 
surface S . In view of the 
Drude-Born-Fedorov constitu- 
tive relations [1], the total 
exterior electromagnetic field 
satisfies in the source-free 

region Q + the modified Max- 
well equations 

VxE‘Jr) = P/ 2 E‘ o (r) + 

v (5) 

+ i«p(|-) 2 H[ 0 (r) 

VxH‘(r) = PKX(r) + 

y, t ( 6 ) 

-iws(^) 2 E^(r) 

where 

/c 2 = w 2 ep , y 2 = i _^ 2p2 (7) 



independently and with dif- 
ferent phase speeds. So, we 
consider the Bohren [9], de- 
composition of E r0 and H r 
into suitable Beltrami fields 
Q, r and Q„ r , as follows 

E r 0 = Ql,t 0 + Q/?,r 0 ( 9 ) 

(10) 

where q = (p/8) 1/2 is the impe- 
dance of the material. 

Because the transformation 
Bohren the relation (8) is 
written : 



' Vx Q t , r ; 




Yl 


0 


1 

lO 

1 


_ V X Q/?,r e _ 




0 


~Yr_ 





with |(3/c| <1, [1] . 

We note that (in contrast 
to the non-chiral case) k is 
not a wave number, but just a 
parameter without physical 
significance which has dimen- 
sions of inverse length and 
a) > 0 is the angular fre- 
quency. The system of rela- 
tions (5), (6) can be written 
as 



or equivalent 

j V X Q t ,r. - Y L %r. 

j^v x Q Rrg = — y R Q R/re 



( 12 ) 



where y L , y R are wave numbers 
for Beltrami fields Q, r , Q„ r 

^■L,r 0 r ^K,r 0 

and are given by, 



Yl = 



k 

1 - Ar(3 



Yr = 



k 

1 + /c(3 



(13) 



~V x E* 
r 0 

VxH' 
r 0 

Beltrami Fields 

As it is well known [1], 
[10], in chiral media LCP and 
RCP waves can both propagate 



respectively. Applying the 

transformation Bohren we 

split the field E[ , H[ into 
the Beltrami fields Q L , Q„ , 

so that the problem (2)-(6) 
becomes equivalent to the 

following problem: We are 



(Y* 



k 2 p iwp 
-iwe /c 2 p 



( 8 ) 
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looking for Q L , Q R in Q + 
such that 

V x Q L (r) = y L Q L (r) , r e Q + (14) 
VxQ R (r) = -y R Q R (r), reQ + (15) 

nxQ t (r) = i^|nxQ s (r), reS (16) 

r x Q t (r) + iQ t (r) = o(i) , r -> oo ( 17 ) 

r x Q R (r)-iQ R (r) = o(i), r ^ oo (18) 



Spherical Beltrami Fields 

So, we consider the Bohren 
decomposition of E^ nc and H“ c 
into suitable incident sphe- 
rical Beltrami fields 

and which have the form 

[6], [7], 

Q^(r|pJ = ^B t (r,r 0 ).p t (19) 

Q^(r|p R ) = ^ R B R (r, r B )-p R ( 20 ) 

where 



B L (r , r e ) = + iw + V x i)h{y L \r - r 0 |) (21) 

B„(r, r 0 ) = ^%y R I + fw - V x I)h(y R |r - r 0 |) (22) 

oTT/ Y R 



and h(x) = e lx /(ix) is the ze- 
roth order spherical Hankel 
function of the first kind 

and I = xx + ijjijj + zz is the 
identity dyadic. We recall 
that B L and B R are the funda- 
mental Green dyadics, for the 
Beltrami fields [1], [11]. 

The constant unit vectors p L 
and p R are assumed to satisfy 
the relations: 

•V P, = f 0 P R = 0 

r 0 xp L = ip L (23) 

r.xp B = - ip R 



The constants A L and A R 
are evaluated so that as the 
location of the point source 
goes to infinity along the 
ray in the direction r 0 , the 
spherical fields degenerate 
into plane LCP and RCP Bel- 
trami fields propagated in a 
direction from r a towards 0 , 
with polarizations p L and 
p R , respectively. Using the 
asymptotic forms: 

|r-r 0 | = r 0 -r 0 -r + 0(i), 

1 0 

p— = - f a + O(-p), (24) 

\r ~ r 0 | r e 

r 0 -> oo 
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in the relations (19) 

(20) we obtain [6], [7], 

Q^(r|p j ) = A J .2e_^K j( -r 0 )- 
'0 

.p. e -^ r + 0(i), (25) 

*i n 



Qf; r :(ripj-r,(^^)p t (3D 

Q^(r|P R ) = ^(£^)-P R (32) 

where 



r e -> w 



for j = L, R , where 

K,(-r o ) = J(I-r 0 r 0 -ir o xl) (26) 

K R (-r 0 ) = |(I-r 0 r 0 +ir 0 xI) (27) 

In view of (23) the dy- 
adics K l and K R satisfy the 
relations 



jp L = K,(-r 0 )-p t 
lp R =K R (-r 0 )-p R 

Hence, if we take 



(28) 



r L =2^(y,I + i vv + Vxl) (33) 

f ' = 4 (/Bi+ ^ VV_Vxi) (34) 



and u = |r-r 0 |. Using asympto- 
tic forms (24) for r oo we 
obtain [6], [7], 



Qj^(r | Pj) = F^(f | PjMYjr) + 

+ o(i) 

where 



(35) 



F^(r I Pj) = 



e~ lyjr ' re 

^(y/ s ) 



K,(r)Pj (36) 



A 

A 



4ny 2 

ky L 

,-iw. inz! 

ky R 



then 



(29) 



lim Q^(r | Pj ) = e 



_ _ 



= Qf (r; -r B , pj) 



(30) 



It is convenient to write 
the incident spherical Bel- 
trami fields as [6], [7], 



are the far-field patterns of 
the point source incident 
Beltrami fields, which sati- 
sfy the relations 

r-F/™(r | Pj) = 0 ( 37 ) 

The scattered electric 
field E* c will be dependent 
on the polarization p L , p R 
and will be have the decom- 
position 



E r s C ( r I Pi./ Pr) = Q4r e ( r I P L , P r ) 
+ Q"r e ( r \PuPr) 



(38) 
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whprp *ZS r I P,' Pr) = h (Yf) + 

(39 ) 

Q"r 8 ( r I P L ' Pr) and *ZS r I P Lf Pr) + O(-i-) 

are the corresponding scat- r 

tered Beltrami fields which The functions g ( f) given by 

have the following behavior, J,e 

when r -» oo [l], [10], [12] 

g ,,r 8 ( r) = ^-K t (f)-J n x (y L V x E“(r') + y 2 E~(r') )e^ f r cfs(r') (40) 

g R , ro (r) = K R (r)- f n x (y R V x E~(r') - y 2 E~(r') )e ^ r 'ds( r') (41) 

are the LCP and RCP far-field where S is the surface of 
patterns, respectively [10], the scatterer £T and n is 

and they are dependent also the outward unit vector on 

on the polarization p L and s . We , also, consider two 

p R . We note that the far- locations for the point sour- 

field patterns g 3r (f) satisfy ae ' a and b from which the 

J,re time-harmonic incident sphe- 



the relations 

fg J , re (r) = 0 
■ f xg L re (f) = -ig L ro (f) 
f xg R re (r) = ig Rir0 (r) 



rical electric waves 



Ef c (r|p„p R ) = 

= Ql”o ( r I Pj + Qr"o ( r I P R ) 



for a = a, b , emanate. E* c and 
the corresponding scattered 
electric waves E* c of form 



Reciprocity and General 

Scattering Theorems electric waves E a of form 

For two vector functions < 38 )' have the Bohren decom - 
u and v we introduce the bi- Position of (9) in terms of 
linear form r 10 1 LCP and RCP Beltrami fields. 

In particular, we have 

{U,V> S = r . 

I £inc _ Qinc _j_ Q 2nc 

= f n (u x V x v - v x V x u)dS- (43) ° 0 ° ( 45 ) 

J v ' ' ' psc qSC qSC 

S ( c o Xl_,o T XR,o 

- 2(3y 2 j n (u x v)dS and the following properties 

s for the fields Q“ c and Q“ 0 

[6], [7]: 
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{Qf"a(r; pj, Qf b (r I Qu %)}s 



2y 2 

k 



J n (Q“ c (r; pj x Q“(r | q L , q R ) )ds{ r) 



(46) 



{QlT (r; pJ , Q"„(r I q R )} s = 0 



(47) 



{Q^a(r;p R ), Q"b(r | q L , q R )} s = - n-(Q^ a c (r;p R ) x Q*'„(r | q L , q*))ds(r) (48) 



{Q R ; a c (r; P R ),QL C b(r-|q„q R )} s = 0 



(49) 



Now a "reciprocity theo- 
rem" for spherical electro- 
magnetic waves in chiral me- 
dia is formulated as follows 
[6], [7]: 

For any two incident 
spherical electric waves 

Ea” c (r I P L , p R ) and Ef c (r|p,,pJ 
of (44) and for any scatterer 
in a homogeneous isotropic 
chiral medium, we have: 



-i V L a 



Qr b ( a iq^q R )-PL + 

ae ^Q'bCalqL, q R )-p R = 



y R 

be- iY ^ 

Yl 

be- ±Y ° b 

V R 



Q“.( b I q u q R )-q L + 

Q R 'a(b I q„ q R ) q R 



(50) 



Let 

Ef c (r|p L ,p R ) and Ef c (r | p L , pj 
be two incident spherical wa- 
ves of the form (44). We de- 
fine spherical far-field pat- 
tern generators for LCP and 
RCP spherical Beltrami fields 
by 



0j, b ( a I q^ q R ) = 



ae 



i Yi a 



1 Vi 



g J .b(a|qb,q R )-K J (f)e 1 ^ a ds(r) - 

S 2 

-Q-: b (a|q„q R )} (51) 



for j = L, R . 

This terminology and defi- 
nition is appropriate because 
when both the observation 
point and the source goes to 
infinity, G Jib (a | q„ q R ) , reduce 

to the far-field patterns for 
an incident plane electric 
wave propagating in the di- 
rection -b and of polariza- 
tions q L for LCP and q R for 
RCP fields. Far-field pattern 
generators in acoustic and 
achiral electromagnetic scat- 
tering have been defined in 
[13] . When the point sources 
are transferred at infinity, 
the generators G., b (a | q L , q R ) are 

transformed into far-field 
patterns scattering spherical 
electrical wave, namely [6], 
[7], 
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limGj b (a|q L ,q R ) = g ■ b (-a|q L ,q R ) (52) 

Using this notation, the 
"general scattering theorem" 
for spherical electric waves 
in chiral media is formulated 
as follows [6] , [7] : 

Let 

E a” c O' I fV P R ) and E“ c (r|p,,pJ 
be two spherical electric wa- 
ves of (44) incident upon a 



scatterer in a chiral medium. 
Then it is valid 

G L,b( a I q L , q R )p l + G R,b( a l q t . q R ) p« + 

+ G L,a(b I P u P R ) q L + G R,a(b I P u Pr) q R = 

= h ^) F-(f)dS(f) + 

Zn s’ Yl 

+ f ^^) FR S ,b(r)t/S(f)} (53) 

? Yr 
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Construction, Measurements and FLOSS 
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Abstract 



The array factor of a both geometrically and electrically 
uniform array is the simple formula for the complex geometric 
progression sum. This fact, although results in the simplest 
of all possible analytical designs, obviously does not in the 
least simplify the complicated practical problem of feeding 
the array elements using multiple driving points. In order to 
begin the examination of uniform linear arrays with a single 
driving point, this paper presents a compact study of the 
end-fed space arrays with application to geometrically uni- 
form, self -standing linear arrays of parallel dipoles. A num- 
ber of test array models were simulated, constructed and 
their radiation pattern was then measured. The experimental 
and computational results were found to be in good agreement. 
The developed software applications are available through the 
Internet as FLOSS Free Libre Open Source Software. 

Keywords 

End-fed, single driving-point, self -standing arrays 



Introduction 

A space array of 1 < k < N 
parallel, arbitrarily shaped, 
dipoles with identical cur- 
rent density distributions, 
has a complex vector radia- 
tion pattern E = AG (by apply- 
ing the "Principle of Radia- 
tion Patterns Multiplication" , 
i.e. by the linear property 
of the volume integral), whe- 
re G is the radiation pattern 
of the first dipole (which 



"generates the array"; the 
"Generator Pattern"), and A 
is the Array Factor, 




i.e. the complex numerical ra- 
diation pattern of N invented 
"isotropic point sources", 
each of current I R and poin- 
ted by the dipole center vec- 
tor R k , with projection R k to 
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the unit direction vector r. 
Linear Arrays have dipole cen- 
ters on a straight line. Ful- 
ly Uniform Linear Arrays are 
both geometrically uniform, 
i.e. the dipoles are equidis- 
tant, and electrically uni- 
form, i.e. the consecutive di- 
pole currents are of equal 
amplitude and constant phase 
difference [1] . 

End -Fed Arrays 

Perhaps, the simplest pra- 
ctical array is the one con- 
structed from a linear two- 
wire transmission line that 
supports the arms of paral- 
lel, linear, symmetrical di- 
poles, vertical to its plane. 
In order to operate the line 
as a balanced one, it is end- 
fed through a coaxial line 
balun from a coaxial conne- 
ctor, which at the same time 
supports the weight of the 
whole, self -standing, array. 
A possible analysis of such 
an array is shown in Fig. 1. 

In this one port linear 
network, the number of cir- 
cuit voltage and current va- 
riables is equal to the sum 
of 2 variables for the input 
port plus 2N variables for 
the dipole ports plus 4(N-1) 
variables for the line seg- 



ment ports, that is a total 
number of 6N-2 variables. The 
number of the linear rela- 
tions between these variables 
is equal to the sum of 3 re- 
lations, for the generator 
dipole 



1 < N 


: V = V = .V 

i ii' 






I = I + . 1 , 

i ii 


(2) 


plus 


2 relations, for 


the 


last 


dipole 




2 < N 


: V = V , 

O N-l N' 






I„ +1=0 

O N-l N 


(3) 


plus 


3 ( N-2 ) relations, 


for 



the intermediate dipoles 



3 < N, 1 < k < N-2 : 
V 



o V k 



k+1 



A+i' 



o^k + ■'‘k+1 + i^k+l ® 



(4) 



plus N relations, for the cou- 
pling between these dipoles 

N 

1 < k < N : V k =Xz kM I, (5) 

M-l 

plus 2 ( N-l ) relations, for 
each one of the lossless 
transmission line segments of 
electrical length (31 and 

characteristic impedance Z 0k 



0 < |31 k t vtt, v = 1, 2, 3, . . . 



ii^k io^k 


_ ^ z 0k 


cospi k 


1 


oi^k oo^kj 


sinpi k 


1 


cospi k 



pl k = vn, v = 2p + 1, p = 0,l,2, . . . : 0 V k = - i V k , 



Ik' 



I, - .1, = 0 

ok l k 



(7) 
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B1 =vtt, v = 2u + 2, u = 0,l,2, . . . : V = V., L + ,1 = 0 (8) 



1 k N 



Tl 

I 






V L 


i l k-l o | 


| i k N-l o 


<r' 







Fig. 1: End-fed linear array 
of linear dipoles 



that is a total of 6N-3 li- 
near relations to solve for 
6N-2 variables, in terms of 
the source voltage V, which 
is considered as a parameter. 
Anyhow, the resulting current 
ratios are clearly indepen- 
dent of V. 

Single Driving -Point 

Self -Standing Linear Array 

The authors' group has li- 
mited available technical re- 
sources for antenna construc- 
tion. This fact restricts the 
experimentation to thin-wire 
self -standing array models of 
a low total weight that is of 
a small total transmission 
ine length and of a small 
number of dipoles. Thus, the 
practical application of the 
analysis was carried-out for 
N = 2, 3 and 4 dipoles only. 
In order to demonstrate the 
procedure in use, the smal- 
lest case of N = 2 dipoles is 



presented in some detail. In 
Fig. 2, the resulting 6N-3 = 9 
linear relations between 6N-2 
= 9 variables + 1 parameter, 
are shown in a compact form, 
for the case of (J1 *vti, v = l, 

2, ... in (6), with each cell 
value to be the coefficient 
of the variable in the first 
row of its column in an im- 
plied summation, while, if 
(31 1 =vn, v = 1, 2, ... then the 

two last rows with the gray 
background have to substituted 
by the rows of Fig. 3 or 4, 
according to the odd or even 
value of v given in (7)-(8). 



Vi 


V 2 


II 


I 2 


dYl 


oVi 


ill 


o 

H 

h* 
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= 
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1 
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l 


1 








-1 
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0 
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1 




-i 


= 


0 
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-1 
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o 
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l 
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-1 




Zn 


Zl2 
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0 




-1 


Z 21 


CM 

r-S 1 
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0 










-1 




$ 

•H 


io z : 




= 


0 












-1 


oi z i 


ooZi 




= 


0 



Fig. 2: The linear system of 
9 relations for arrays with 
N = 2 and pi x t vtt, v = 1, 2, ... 

Obviously, the complexity 
of the expressions increase 
with the number of dipoles, 
from the simple, of 2 dipoles 
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I 




-12oo Z l io Z l Z 22 



io Z l"*"io Z l Z 12 oo^-lV 



( 9 ) 



in which the equality of self 
and mutual impedances have 
been taken into account, re- 
sulting from the system of 
Fig. 2, or of the simplest 
I 2 /I 1 = +l or I 2 /I 1 = -l of Fig. 

3 or Fig. 4 respectively, to 
the most complex one for 4 
dipoles, shown in Fig. 5, 
which covered about one and a 
half A4 page. 











1 


1 








= 


0 














-1 


1 




= 


0 



Fig. 3: Last rows replacement 
for v = 2p + 1, p = 0, 1, 2, ... 











-1 


1 








= 


0 














1 


1 




= 


0 



Fig. 4: Last rows replacement 
for v = 2p + 2, p = 0, 1, 2, ... 




i^l o^l i.^2 o^2 i^3 0^3 



Fig. 5: Equivalent circuit of the 4 dipoles linear array 



Three Visual Fortran ap- 
plications were developed for 
the computation of current 
ratios. The GUI application 
form for N = 4 dipoles, is 
shown in Fig. 6 [2], [3]. In 
this form, the input data are 
the N-l distances between di- 
poles, the dipole radius and 
length, the length, the char- 
acteristic impedance and the 



velocity factor of each tran- 
smission line segment. In ad- 
dition to current ratios, 
each application produces the 
text files needed by the 
[RadPat4W] application of a 
FLOSS mini-Suite of tools, 
which plots in Virtual Real- 
ity the Generator, Array Fac- 
tor, and Array 3D radiation 
pattern, as well as, their 2D 
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main-plane cuts [4]. The de- 
duced formulas for the deter- 
mination of the 1 + 2 + 3 = 6 
current ratios (I k /I 1 ), were 

mechanically verified using 
Mathematica . 

Array Design, Construction 

and Measurements 

A number of dipole array 
designs were carried out us- 
ing the developed applica- 
tions and the [RichWire] si- 
mulation program, which is a 
fully analyzed, corrected and 
redeveloped edition of the 
original Moment Method thin- 
wire computer program [4], 

[5] . Eleven arrays were fi- 
nally constructed and tested. 
The results for the current 
ratios of four selected ar- 
rays are shown in Tab. 1. The 
arrays were designed for op- 
eration at the frequency of 
1.111 GHz. A two-wire trans- 
mission line, of Z 0 = 200 £2 and 

velocity factor vf = 1 was con- 
structed to feed the dipoles 

[6] . This balanced line was 
then connected to an unbal- 
anced 50 £2 type-N/F base con- 
nector through a 4:1 balun 
made from a segment of RG- 
174U coaxial cable (Z 0 = 50 £2, 

vf = 0.66) with total length 
A/2. 

The arrays were constru- 
cted by bare copper wire of 1 
mm (0.0037A) radius and they 
are self -standing using an 
orthogonal piece of a two- 



sided printed board (3 cm x 
4.48 cm), on which the two- 
conductor line was soldered. 
A few Teflon spacers of low 
relative dielectric constant 
(= 2) were fabricated to fix 
the distance between the two 
line wires at 5.5 mm. The 
measurement system consists 
of a 50 £2 Vector Network Ana- 
lyzer external to an anechoic 
chamber [7] . Each array under 
test was azimuthally rotated 
around its three main axes, 
by a 360° built positioner, 
under the developed software 
control of a built hardware 
controller. The stationary an- 
tenna was a UHF standard gain 
antenna [8]. Fig. 7 and Fig. 
8 show the results for the 4 
test arrays, in 3 groups of 
rows, as follows: 1st group: 
The screen captures of the 
produced Virtual Reality ra- 
diation patterns in dB for 
(a) Generator, (b) Array Fac- 
tor, and (c) Dipole Array - 
2nd group: The constructed 

model, and 3rd group: The 2D 
radiation pattern cuts by (d) 
xOy, (e) yOz, and (f) zOx 
main-planes. In Fig. 7, the 
measurements for the array of 
N = 2, dipoles with equidis- 
tance d = 0.325A, were carried 
out with and without balun. 
In Fig. 8, the array of N = 2 
dipoles with equidistance d = 
0.633A, has been designed to 
exhibit the maximum radiation 
pattern direction off the 3 
main-axes . 
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Fig. 6: GUI for the analysis of a linear array 
of N = 4 dipoles 



Tab. 1: The experimental array characteristics and 

current ratios 



Dipoles 


2 


3 


4 


2 


L/A 


0.5185 


0.5 


0.5 


0.44 


d/A 


0.325 


0.325 


0.325 


0.633 


Ratios 


V 1 ! 


V 1 ! 


V 1 ! 


V 1 ! 


V 1 ! 


V 1 ! 


V 1 ! 


IVM 


0.534 


0.205 


0.286 


0.280 


0.115 


0.226 


0.527 


^I k /Il[°] 


-80.5 


-114.9 


173 


-115.9 


150.4 


68.1 


49.5 
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Conclusion 

Although, the unavoidable 
mechanical supporting ele- 
ments which exist in the ane- 
choic chamber in the neigh- 
borhood of the antenna under 
test may affect the radiation 
pattern measurements, the ob- 
served differences in the 4 
of the 12 patterns between 
the analysis and simulation 
results, on the one hand, and 
the measurements, on the 
other, have to be charged re- 
spectively: (1) On a cone-cut 
of the array radiation pat- 
tern, instead of the expected 
yOz main-plane cut, in Fig. 
7(e), N = 2 and Fig. 8(e), N = 
4, (2) On an inclination of 
the rotation axis relative to 
the expected linear polariza- 
tion measurement plane, in 
Fig. 8(d), N = 2, and (3) On a 



loosed connection during the 
array rotation, in Fig. 7(d), 
N = 3. These conclusions are 
amplified by the careful stu- 
dy of the corresponding Vir- 
tual Reality space radiation 
patterns. Therefore, under the 
given measurement circumstan- 
ces, the experimental and com- 
putational results were found 
to be in good agreement and 
no attempt was made to modify 
any design or repeat any mea- 
surement . 

The results for the single 
driving-point self -standing 
fully Uniform Linear Arrays, 
i.e. those including the elec- 
trical uniformity, along with 
their application to the con- 
strained pattern design will 
be presented in a future pa- 
per (Part 2) . 
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Abstract 

The aim of this paper is the study of various parameters 
affecting the ageing of transformer oil in distribution 
transformers of 20/0.4 kV. Fifty (50) samples of oil were 
taken from such transformers. The transformers function in 
the major Athens area, Greece. Parameters, such as breakdown 
strength, oil color, humidity, interfacial tension and tan5 
were taken into account. Transformer ageing and lifetime are 
strictly related to the rate of ageing of the whole insulat- 
ing system, and mainly of the oil. 

Absence of transpiring system with silica gel has as a 
result the increase of oil humidity. The lengthy use of 
transformers under heavy load, and consequently under high 
temperatures, is a main factor for oil ageing and oxidation. 
Especially in the last few years, the increased loads re- 
quired by the transformers which have to do also with the 
climatic changes, have as a result the additional stressing 
of the oil. Various arcs, resulting from short circuits in 
the network, have as a consequence the production of gases 
and sludge. Such gases and sludge influence in a negative 
way the insulating properties of the oil as well as its rate 
of ageing. The role of silica gel is stressed. 

Keywords 

Transformer oil, distribution transformer, diagnostic me- 
thods, breakdown voltage, dielectric strength 
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Introduction 

Transformer oil is a most 
important component in trans- 
formers [1] . Transformer oil 
ageing has been studied as 
well as the parameters and 
factors influencing its be- 
havior under a variety of 
stresses (electrical, ther- 
mal, etc.) [2] - [5] . Faults in 
power and distribution trans- 
formers are rare (of the or- 
der of 1% - 2% per year), but 
when they occur, they have 
very serious technical and 
economic consequences. They 
may even lead to dangerous 
situations for human life and 
the environment. Main factors 
affecting the acceleration of 
ageing of the insulation of a 
transformer are humidity, high 
temperature, oxidation and the 
acidity of its oil. The role 
of chemical byproducts of the 
transformer insulation is al- 
so very important for its 
ageing. There is no single 
measurement able to give a 
whole picture of the state of 
a transformer. It must be no- 
ted that the variety of the 
diagnostic methods used, does 
not have as an aim the pre- 
diction of the useful remain- 
ing lifetime of the insula- 
tion, but it tries to reveal 
the increasing probability of 
faults and the corresponding 
decrease of the insulation 
reliability. Needless to say, 
that reduced reliability im- 
plies also reduced remaining 



insulation lifetime. 

With this in mind, in the 
context of the present paper, 
diagnostic methods were used 
in order to study the state 
of distribution transformers 
of the major Athens area. The 
transformers investigated we- 
re 20/0.4 kV. The whole work 
was carried with the aid of 
Public Power Company (PPC) of 
Greece, and more specifically 
with the aid of PPC Trans- 
former Division. Oil was ta- 
ken from 50 distribution 
transformers . 

Diagnostic Methods 

Warning signs about the 
state of a transformer are, 
among others, a big increase 
of partial discharges (» 
2500 pC), a visible deterio- 
ration because of foreign me- 
tallic and carbon particles, 
the presence of humidity in 
the solid insulation about 3- 
4 % and the presence of 
sludge . 

Several diagnostic methods 
were used in order to see the 
quality of the transformers 
in question. The characteri- 
zation of the oil color (DIN 
51517 - ASTM D155 ) was per- 
formed through a device 
(chromometer) including stan- 
dard glass disks and two 
glass jars with lid. The con- 
trol of dielectric strength 
was measured by a typical 
Foster test cell, according 
to IEC 156/95 (Fig. 1). 
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The control of humidity in 
the oil was measured by a 
Metrohm - 684 KF Coulometer, 
which consisted of a glass 
container with a stirrer ti- 
tration in which the reagent 
from container storage is ad- 
ded. The device is fully au- 
tomated and once the experi- 
menter gives the settings, it 



measures the moisture content 
of the oil. The measurements 
were performed according to 
IEC 814. The control of in- 
terfacial tension ( ASTM D971 
- 91) was performed via a 
tensimeter, which gives the 
value in dynes per centimeter 
in a direct reading (Fig. 2). 




Fig. 1: Device for the measurement of dielectric strength 
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The device that performed 
measurements of tan5 and of 
resistivity, is the BAUR-DTL 
fully automated device for 
measuring dielectric losses 
of oils. Such a system has a 
fully automated process for 
measuring dielectric loss, 
relative dielectric constant 
and resistivity. The measure- 
ments were performed with a 



system counting tan5 values 
with maximum accuracy from 
0.00001 to 4.0, measured ac- 
cording to IEC 247. 

It is true that no single 
diagnostic method can give 
full information as to the 
state of a transformer. The 
aforementioned methods may 
give a better picture of its 
state . 




Fig. 2: Tensimeter Cenco du Nouy 
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Experimental Results 

The sampling was performed 
with due care, and according 
to general standard practice, 
i.e. during sampling there 
must not be any dust and hu- 
midity in the nearby space, 
the sampling cells must be 
clean and dry and they must 
be washed with oil from the 
transformer which is to be 
checked, the samples must be 
protected from light and they 
must be taken while the oil 
is hot but due time should 
pass before the various par- 
ticles settle down. Sampling 
was done from the bottom of 
the transformers. The glass 
cells used for oil sampling 
were big enough so that they 
contained enough oil volume 
for additional measurements, 
if needed. No measurement was 
performed based on only one 
value or sample. All results 
must be verified with addi- 
tional samples and measure- 
ments. Color measurements of 
the oil samples were based on 
specification DIN 51517 
ASTM D155 (Fig. 3). Dielec- 
tric strength measurements 
were based on specification 
IEC 156/95. Humidity measure- 
ments were based on specifi- 
cation IEC 814. Interfacial 
tension measurements were ba- 
sed on specification ASTM 
D971-91 . Loss factor measure- 
ments were based on specifi- 
cation IEC 247. 

Transformer oil is charac- 



terized as good, if its color 
is not dark, its breakdown 
voltage more than 40 kV, its 
humidity less than 10 ppm, 
its resistivity more than 3 
GD-m, its loss factor 0.1, 
its interfacial tension more 
than 28 mN/m. It is accept- 
able if its breakdown voltage 
is between 30 to 40 kV, its 
humidity between 10 and 25 
ppm, its resistivity between 
0.2 and 3 GQ-m, its loss fac- 
tor between 0.1 and 0.5, its 
interfacial tension between 
22 and 28 mN/m. In the case 
of an acceptable oil, sam- 
plings should be carried out 
more frequently, and the con- 
tent in foreign particles and 
contained water should be 
controlled. If the breakdown 
voltage and humidity are near 
the limit values, then the 
oil should be filtered and 
cleaned. An oil is poor when 
its color is dark (the dark 
color is an indication of 
pollution or ageing), its 
breakdown voltage less than 
30 kV (in this case the oil 
should be replaced or thor- 
oughly cleaned), its humidity 
more than 25 ppm (in such a 
case the oil should be re- 
placed or thoroughly clea- 
ned), its resistivity less 
than 0.2 GQ-m, its loss fac- 
tor more than 0.5, its inter- 
facial tension less than 22 
mN/m (in this case a check 
should be carried out for the 
presence of sludge). 
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Generally it is suggested 
that in case one of the pa- 
rameters mentioned (color, di- 
electric strength, humidity, 
interfacial tension, loss fa- 
ctor, resistivity) is accept- 
able, a more frequent sam- 
pling should be done in order 
to better monitor the oil 
quality. In the case any of 
the mentioned parameters is 
poor, cleaning or replacement 
of the oil is suggested. In 
any case, it is good to con- 
trol all parameters concerned 
in order to have a clear pic- 
ture of the state of the oil. 

From the 50 oil samples 
investigated (taken from 50 
distribution transformers), 3 
of them showed breakdown vol- 
tage lower than 30 kV (i.e. 
percentage of 6%), 45 samples 
indicated humidity higher 



than 25 ppm (i.e. percentage 
90%), 28 samples presented 
interfacial tension lower 
than 15 dynes/cm (percen- 
tage 56%), whereas 44 samples 
showed a color higher than 1 
1/2 (percentage 88%) (Fig. 
4) . It was evident from the 
data - mentioned in detail in 
[6] - that the more the number 
of years in service, the more 
the number of samples with 
lower dielectric strength. 
This is due to the fact that 
with the years there is also 
an increase of humidity in 
oil as well as an increase of 
the oxidation byproducts. 
Trans formers which are in 
cities or in urban localities 
do not usually suffer light- 
ning strokes, consequently 
the only stresses come from 
switching overvoltages and/or 
some high currents. 




Fig. 3: Chromometer Hellige Comparator 
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Fig. 4: Typical glass colored disc 



As the years pass by, hu- 
midity increases because fluc- 
tuations in temperature cause 
an increase of humidity en- 
tering the main oil volume. 
In distribution transformers 
having silica gel, the latter 
looses its absorbing property 
with an increasing number of 
years in service. Consequently 
silica gel should be replaced 
at regular intervals. The in- 
creasing number of years in 
service influences also the 
oil oxidation. Oxidation by- 
products in turn affect tan5. 
The loss factor, however, is 
not influenced from the exis- 
tence of humidity in the sam- 
ples and for this reason is 
not a criterion for its exis- 



tence. Interfacial tension is 
reduced with the years in 
service. This is due to the 
increasing quantity of humid- 
ity as well as to the bypro- 
ducts of oxidation. There is, 
however, a number of samples 
with rather reduced interfa- 
cial tension although the 
number of years in service 
was not that big. This may be 
due to the quality of the oil 
used or to the fact that it 
contained only a small quan- 
tity of anti-oxidants. Fi- 
nally, the oil color changes 
with the years in service. A 
change of color may also be 
due to the overcharging of a 
transformer . 
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Discussion 

In this work, sampling from 
50 transformers of 20/0.4 kV 
was carried out. Sampling and 
measurements were performed 
as suggested by the various 
international specifications 
[7] . All transformers were in 
the major Athens area. It 
should be noted that the trans- 
formers were not all from the 
same manufacturer. There were 
transformers from various ma- 
nufacturers, both from Gree- 
ce and abroad. This is one of 
the complications of the Greek 
Electricity System and cer- 
tainly, this fact hinders any 
thought for a statistical ap- 
proach. In case of poor oil 
color, this may be a symptom 
of humidity. Poor oil dielec- 
tric strength may imply that 
the oil should be cleaned or 
replaced. If the humidity le- 
vel is high, the oil has to 
be again either cleaned or 
replaced. If the interfacial 
tension of oil is acceptable, 
this means that a control for 
the possible presence of by- 
products may be performed. 
Acceptable values of oil re- 
sistivity may be tolerated. 
In most cases, tan5 was good. 
The parameters which changed 
most were the color, the oil 
dielectric strength, the lev- 
els of humidity and the oil 
interfacial tension. From the 
above, it is evident that one 
cannot pronounce any verdict 



on the quality of a trans- 
former oil based on only one 
parameter . 

It is evident from the re- 
sults that there was absence 
of a transpiring system with 
silica gel in most investi- 
gated transformers. This ab- 
sence has as a result that 
the humidity of the atmos- 
phere is in contact with the 
oil. The functioning of a 
transformer under heavy load, 
and consequently its func- 
tioning at high temperatures, 
results in an acceleration of 
its ageing. High temperatures 
cause oxidation of the oil. 
In addition to that, in the 
last years, increased loads - 
as a consequence of climatic 
changes (such as increased 
temperature and humidity) - 
are required and they result 
in the stressing of the oil. 
Arcs, because of short cir- 
cuits, create gaseous byprod- 
ucts and sludge. The latter 
influence in a negative way 
the dielectric strength of 
the oil. Another point which 
should be considered is that, 
with the continuing effort 
from the side of the manufac- 
turers for the reduction in 
the size of the transformers, 
the oil tends to be more eas- 
ily thermally stressed. This 
results in an accelerated 
rate of rise in temperature 
and an accelerated rate in 
its ageing. 
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Generally speaking, the 
frequent oil sampling helps 
in pointing out potential 
trouble spots in the network. 
It goes without say that, as 
in this paper, a variety of 
measurements may give a bet- 
ter picture of the actual 
state of the transformers ra- 
ther than isolated measure- 
ments or measurements confi- 
ned only on a single quantity 
(e.g. dielectric strength). 
Ageing is a complex process 
and, as is well known, a va- 
riety of factors play a role. 
Ageing factors contribute not 
only to the ageing of the 
oil, but they may also inter- 
act among themselves, render- 
ing thus the prediction of 
lifetime even more difficult 
[8], [9]. Breakdown of trans- 
former oil is also a complex 
phenomenon determined mostly 
by the operating conditions 
[ 10 ], 

As a way out of the vari- 
ous problems with transformer 
ageing, the construction of 
additional substations is pro- 
posed. This will reduce the 
load of each transformer. Need- 
less to say that, the contin- 
uous checking of the whole 
distribution network and the 
requirement from the manufac- 
turers to observe all the in- 
ternational norms and speci- 
fications, are necessary con- 
ditions for the good func- 
tioning of the whole system. 



Special attention should 
be paid to the way one gets 
the oil samples. Sampling 
should be done with due care 
and care should be taken in 
order not for the sample to 
contain water. The samples 
should be taken from the 
transformers when the oil is 
still warm. Foreign particles 
should be allowed to settle 
at the bottom, so that no un- 
due influence can be recorded 
in the various measurements. 

In the case the results 
are contradictory, one has to 
repeat the sampling and the 
corresponding measurements. 
It is true that the manufac- 
turers, try to reduce the pro- 
duction cost of the trans- 
formers and their respective 
materials by reducing the si- 
ze of the transformers. This 
on one hand has as a result 
the possibly higher break-down 
strength of the oil, on the 
other hand, however, has as a 
result the faster increase of 
the oil temperature, and con- 
sequently the faster ageing 
and stressing of the mechani- 
cal parts of the transformer . 
Anyhow, a better maintenance 
of the distribution transfor- 
mers should be followed as 
well as the placing of more 
transformers of larger power 
in localities with an in- 
creasing load demand. 

The discussion as well as 
the conclusions of the pre- 
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sent paper are in line with 
previously published work 

[11] , [12]. Such work was 
carried out on transformers 
of 20/0.4 kV as well as of 
150/20 kV. A statistical ap- 
proach cannot be offered for 
the time being, as it was 
done before [13], for the 
reasons given above. 

It is to be hoped that the 
next steps in our work will 
be a classification of the 
distribution transformers ac- 
cording to manufacturer, ac- 
cording to their power (kVA) 
and according to the previous 
history (which includes years 
in service, all faults, short 
-circuits, lightning strokes, 
switching operations etc.). 
Thereafter, an effort will be 
made in order to statisti- 
cally analyze the data. 

Concluding, it is fitting 
to note that the present pa- 
per does not claim any origi- 
nality regarding any new di- 
agnostic techniques. This pa- 
per is an application of al- 
ready existing and well known 
specifications. It is a con- 
tinuation of an effort star- 
ted few years ago, as men- 
tioned in [7], [11], [12]. In 
this respect, the present pa- 
per presents new results, al- 
though the methods used are 
the same as in [6], [7], [11], 

[12] . In other words, this 
paper presents the state of 
transformer oil in some dis- 



tribution transformers of a 
given area at a given moment 
in time. Needless to say that 
continuous monitoring is 
necessary in order to check 
and confirm the transformer 
oil quality. Continuous moni- 
toring, however, is often not 
enough. If possible, it should 
go together with fault dete- 
ction, namely partial dischar- 
ge (PD) detection. PD dete- 
ction can help preventing a 
defect from developing into a 
fault, e.g. a short circuit 
[14], [15], [16]. Nevertheless, 
condition monitoring as pre- 
sented here, should be done 
on a regular basis. 

Conclusion 

Transformers of 20/0.4 kV 
from the major Athens area 
were investigated. Various pa- 
rameters (oil color, dielec- 
tric strength, humidity, in- 
terfacial tension, tan5 and 
resistivity) of the oil were 
studied with standard methods 
based on international speci- 
fications. No verdict on trans- 
former oil quality can be ba- 
sed on only one parameter. A 
multitude of parameters is 
needed in order to pronounce 
a correct verdict on oil qua- 
lity. The continuous control 
and monitoring of the distri- 
bution transformers is neces- 
sary in order to avoid prob- 
lematic situations. The role 
of silica gel is emphasized. 
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Silica gel should be used in 
all transformers, so that hu- 
midity can be absorbed. The 
quantity of humidity increases 



with the number of years in 
service, consequently, distri- 
bution transformers should be 
checked at regular intervals. 
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The thirteen-leaf is a symbol for a life tree leaf. 
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"Warrior with a phiale", ca. 480-460 BC, 

by Berliner Maler, 
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